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Abstract. We consider a simple model for the unfolding of RNA tertiary structure under dynamic loading. 
The opening of such a structure is regarded as a two step process, each corresponding to the overcoming 
of a single energy barrier. The resulting two-barrier energy landscape accounts for the dependence of 
the unfolding kinetics on the pulling rate. Furthermore at intermediate force, the two barriers cannot be 
distinguished by the analysis of the opening kinetic, which turns out to be dominated by a single macro- 
barrier, whose properties depend non-trivially on the two single barriers. Our results suggest that in pulling 
experiments on RNA molecule containing tertiary structures, the details of the single kinetic barriers can 
only be obtained using a low pulling rate value, or in the high force regime. 

PACS. 87.14.Gg DNA, RNA - 82.37.-j Single molecule kinetics 



1 Introduction 

The study of structural properties of biological molecules 
has received a boost by the introduction of techniques 
allowing for the manipulation of single molecules. For ex- 
ample, the study of folding and unfolding of nucleic acids 
can now be performed by applying a controlled force on 
the free end of a single strain of a molecule. In this situa- 
tion, the opening of the Watson-Crick pairs leads to what 
has been called the unzipping of the molecule. Unzipping 
in DNA is similar to some steps involved in DNA repli- 
cation and in its translation into mRNA, and has been 
the subject of several theoretical and experimental inves- 
tigations [HE]. On the other hand, RNA unzipping ex- 
hibits further complications, since a given single stranded 
RNA molecule can exhibit a complex secondary structure 
(matching pattern between complementary bases) and ter- 
tiary structure (three-dimensional conformation). Thus, 
single-molecule unzipping experiment can yield informa- 
tion on the secondary structure of RNA molecules |3]0J 
I^H^HTj. In ref. [S] in particular, the role of the secondary 
structure intermediates in the folding/unfolding experi- 
ments is discussed, and it is shown that such intermediates 
can be responsible for the slowing down of the kinetics. 
Moreover, the response of complex RNA structures to ap- 
plied mechanical forces can be analogous to the responses 
of RNA during translation or export from the nucleus 0]. 
However, in such cases, the tertiary RNA structure plays 
an important role. It is commonly believed that the break- 
ing of RNA tertiary structures or their formation are the 
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time limiting processes in the molecule unfolding or refold- 
ing, respectively [3]|H]- In the unzipping experiments, ter- 
tiary contacts may lead to the appearance of kinetic bar- 
riers where the unzipping process momentarily stops 
|I] • In this situation, overcoming these barriers bears some 
analogy with the breakup of molecular adhesion studied in 
a number of experiments 9 ,1Q| which have been recently 
reanalyzed theoretically [TT1ll2lll3j . In the present paper, 
we wish to understand some features of the experiments 
on RNA unzipping in the light of this theory, in order to 
highlight what information one is able to collect from the 
unzipping kinetics on the position and the height of the 
barriers due to the tertiary structure. Although in the case 
of the folding of large complex RNA molecules, the forma- 
tion of kinctically trapped intermediates in the secondary 
structure can play an important role in the slowing down 
of the process |14II15II16| . here we will specifically consider 
the role of the tertiary structures in the unzipping ex- 
periments. Such role was stressed in refs. [30], where the 
mechanical unfolding of RNA molecules were performed. 
In those experiments the tertiary contacts could be re- 
moved by changing the solution the RNA was immersed 
in. The removal of the tertiary structures corresponded 
to the disappearing of the kinetic barriers, and further- 
more the folding/unfolding processes became reversible. 
This kind of experiments indicate that the tertiary struc- 
tures are responsible for the kinetic arrest of the structural 
rearrangement of RNA molecules under tension, although 
the tertiary structures are much more brittle compared to 
the secondary ones, i.e they break in consequence of small 
deformations. 
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2 The model 

A typical unfolding experiment consists in holding one 
molecule's free end in an optical trap while the other free 
end is pulled at constant velocity 00]. This induces a 
pulling force on the molecule that increases linearly with 
time f(t) = rt, where r is the loading rate. The typical 
output of such experiment is a force-extension curve where 
the monotonic increase of the force is interrupted by a 
number of plateaus revealing the unfolding of portions of 
the molecule. 

The unfolding of the molecule is a stochastic process, 
which depends on the pulling force rate r, the actual value 
of the pulling force and the microscopic details of the 
molecule. Therefore the unfolding of an RNA molecule, 
i.e., the successive breaking of its molecular links, can 
be viewed as a succession of thermally activated escape 
processes over a set of energy barriers, each representing 
one or more molecular links 8,14,15 . Within this pic- 
ture, the resulting energy landscape can be considered 
one-dimensional, since the experimental set-up singles out 
a well-defined direction, which is the pulling direction. 
The landscape energy barriers are thus located at increas- 
ing distance along the pulling direction. This experiment 
can be repeated for several values of r, in such a way 
that the different force-extension curves can be mapped 
onto a single curve of the breaking force as a function of 
the loading rate. From this curve, much information on 
the microscopic details of the molecular bonds can be ob- 
tained ptTU] : in particular the typical length and energy 
of such bonds. Another quantity, which can be obtained 
by the force-extension characteristics, is the fraction 4> of 
molecules which remain folded as a function of the time 
or the force, and which also depends on the parameter 
r. Sampling the fraction of folded molecules at different 
times provides a more direct insight into the kinetic barri- 
ers which slow down the unfolding process. Using Evans' 
results for the case of pulling experiments on a single 
molecular link jl()| . which can be represented by a sin- 
gle kinetic barriers of height AE and position Ax along 
the reaction coordinate, the fraction <fr(t) of bound links 
at time t is given by 



4>(t) = exp 



at cone H y > 
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where luq is the attempt frequency, which depends on the 
microscopic details of molecular linkage, and (3 = l/ksT. 

In the large force regime (/ > 10 pN), one expects that 
the quantity q(f) defined as 



q(f) = in 
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will be a linear function of /: 



q(f) ~ In 
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and will exhibit no dependence on r. 

This expression has been exploited by Liphardt et al. |3] 
to characterize the tertiary structure of an RNA molecule: 
by mechanically pulling on the P5abc domain of the Tetrahy- 
mena thermophila ribozyme, the fraction of folded molecules 
as a function of the force (time) has been determined. The 
authors then make the hypothesis that the tertiary struc- 
ture can be described as a single kinetic barrier, which 
hinders the molecule from unzipping. Using eq. they 
give an estimate for the two characteristic parameters of 
the barrier, i.e. the zero-force transition rate, defined as 



fc = cj exp (— fiAE) 



(4) 



and the barrier position along the reaction coordinate Ax 
(in that case and in the following in this paper, the elonga- 
tion of the molecule needed to break the bond). However, 
the data for q(f) showed in that paper, exhibit a depen- 
dence on r even if not distinct, in contrast with eqs. (j2J 
and ©• 

In the present paper we will argue that at least an- 
other kinetic barrier has to be considered in the energy 
landscape of the tertiary structure in order to account 
for the dependence of q(f) on r. Furthermore, in a re- 
cent paper, Bartolo et al. ^3] showed that in the more 
complicated case of several molecular links, which can be 
represented by a set of N kinetic barriers along a one- 
dimensional unbinding path, the unbinding force plotted 
as a function of logarithm of the pulling rate (In r) appears 
as a succession of straight lines whose slopes are given by 
the distances between the adjacent maxima and the min- 
ima of the energy landscape. Thus, we expect that for the 
simple two-barrier energy landscape here considered, the 
plot of q as a function of / = rt will exhibit more than 
one straight line, each corresponding to a different escape 
route from the folded to the unfolded state. 



E 




Fig. 1. Schematic plot of one- dimensional energy landscape 
with two energy barriers. The full line corresponds to the land- 
scape with no external force applied E(x), while the dotted 
curve corresponds to the time-dependent landscape E(x) — rtx 
at t > tA (see text and eq. ® for the definition of Ia)- 



Basically we assume that the tertiary structures of a 
\ ' RNA molecule can be described by more than one energy 
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barrier with different width and height: the simplest case 
is the one depicted in fig.H where a landscape E(x) char- 
acterized by two barriers is responsible for the slowing of 
the unfolding process, and the state x = corresponds 
to the unperturbed structure. The opening of the tertiary 
structure can thus be regarded as a two-step process: when 
the molecule extension is of the order of xa, as a conse- 
quence of thermal denaturation or pulling force, a first set 
of tertiary interactions are broken, and the system jumps 
into the local energy minimum at x a . Later, as the length 
xb is reached, the remaining tertiary interactions break 
down. The choice of a two-barrier landscape has been also 
suggested to us by the following experimental observation: 
in RNA pulling experiments, a small but finite fraction of 
molecules under dynamic loading unfolds with two suc- 
cessive rips jnj, revealing the existence of an intermedi- 
ate state between the completely folded and the unfolded 
ones. 

When force is applied, the energy landscape changes as 
E(x, t) = E(x) — rtx, so that the outermost barrier height 
decreases faster than the inner one: the unfolding kinetics 
is thus characterized by two crossover times which will 
be respectively indicated with tA and t a in the following. 
For times longer than tA, the inner barrier becomes the 
dominant one and is mainly responsible for the slowing 
down of the unfolding process with respect to the / 
( j case. Eventually, at times longer than t a (with t a > 
tA), the well (x a , E a ) disappears and the unfolding simply 
results from the overcoming of the first barrier (xa, Ea)- If 
we make the further assumption that the energy barriers 
and wells sketched in fig. ^ are sharp, so that the barrier 
and well positions remain essentially constant with time, 
the crossover time tA is given by 



t A 



E A 



and t a is given by 



t(xb - xa) 



Eb — E a 

r(x B - X a ) 



(5) 



(6) 



The last equation also defines a crossover force f a — r -t a . 
Using a Kramer formalism, and the notations shown in 
fig. ^ we can write the instantaneous rate for the transi- 
tion from one of the two minima (x = 0, x a ) over the two 
corresponding energy barriers (x — xa,xb)'- 

k Q ^ a = luq exp [-0 (E A - f(t)x A )\ , (7) 
k a ^a = w exp [-(3(E A - E a - f(t)(x A - x a ))\ , (8) 
K^oo = uj q exp [-0 (E B - E a - f(t)(x B - x a ))} , (9) 

where x = oo indicates the completely unfolded state. We 
assume furthermore /coo^o.a = 0, i.e., once unfolded the 
system never folds back, as observed experimentally [3]. 
Let Po(t) and p a (t) be the probabilities that the system 
is in the state or a, respectively. The time evolution of 
this quantities is described by the following differential 



equation system 

Po = -ko^aPa + 0(t a - t)k a ^oPa , (10) 

Pa = 0(t a - t)k Q ^ a P0 - 0(t a - t)(k a ^ + k a ^oo)Pa 

-e(t-t a )u , (ii) 

where 9(t) is the Heaviside step function. The step func- 
tions 9 have been included in eqs. (jlOJl and 1)11(1 in order 
to take into account the disappearing of the well (x a ,E a ) 
from the system energy landscape at t — t a - It is worth to 
note that at such crossover time, a significant fraction of 
molecules might be accumulated in the x — Xa state, as 
a result of the system evolution at previous times. Thus 
we assume that at time t > t a , the escape rate of the 
molecules which are still in the state x = x a , is deter- 
mined by the molecular attempt frequency loq alone. For 
a given set of characteristic parameter 

S = {LL>o,x A ,EA,x a ,E a ,XB,E B } , 

given the initial values Po(0) = 1 and p a (0) = 0, such 
a system can be solved numerically. As mentioned above, 
the molecule never folds back, once it has been completely 
unfolded, therefore the quantity 



<p(t)= Po (t)+p a (t) 



(12) 



defines the probability that the system is still folded at 
time t, either completely in the x = state or partially in 
the x = a state. 

In order to obtain a reliable estimate of the param- 
eter set S, we consider the results of the above cited 
experiment: applying eq. © to the openings of a sim- 
ple RNA molecule tertiary structure, a zero-force tran- 
sition rate fco — 2 x 10 _4 s _1 and a difference between 
the unfolded tertiary structure length and its transition 
(breaking) length Ax ~ 1.6 nm have been obtained 
The zero-force transition rate ko = 2 x 10 _4 s _1 is re- 
lated to the attempt frequency ujq and to the overall en- 
ergy barrier AE of the single barrier picture via eq. 
In order for the tertiary structure to be the dominant 
impedance to the molecule unfolding, the involved energy 
barriers have to be greater than the well known RNA base 
pair energies, which are of the order of a few k-gT. Thus, 
if we suppose that the energy barrier AE is of the or- 
der AE ~ 10 k-oT, from the above cited result and from 
eq. |@}, we obtain ujq ~ 4.4 s _1 . Our estimate for the 
tertiary structure energy is in agreement with the values 
shown in |17j . where combining numerical computations 
with experimental techniques, the energy of tertiary inter- 
actions in a simple RNA molecules was found to range be- 
tween 6 k^T and 13 k^T . On the other hand, a direct mea- 
surement of the attempt frequency ojq in folding/unfolding 
experiments is rather difficult. An indirect estimate can be 
obtained by pulling the RNA molecule at constant force, 
when the force value is within an interval of a few pNs 
around the unfolding force. The molecule then hops back 
and forth between the folded and unfolded state with a 
frequency between 0.05 s _1 and 20s -1 [3], which depends 
on the actual value of the force. It can be assumed that, at 
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the unfolding force, the effective energy barrier vanishes 
and the hopping rate yields a rough estimation of the mi- 
croscopic attempt rate. Thus our estimate for the attempt 
frequency u>q is in agreement with those experimental re- 
sults. 

Here we take xa = 0.6 nm, x a = 0.8 nm, and xb = 1.8 
nm: compared to the above cited result, our choice corre- 
sponds to two successive openings which occur at elonga- 
tions whose sum is equal to the single step picture elon- 
gation Ax = 1.6 nm. The set of parameter values that we 
will consider in the following is thus given by 

S = {w =4.4s-\x A =0.6nm,£U=10fc B r, (13) 
a; a =0.8 nm, E a —6 k^T, xb—1-8 nm, Eb=16 &bT} , 

with T = 300 K. Note that with this choice for the energy 
barrier height, after the first opening, the system has to 
overcome another barrier of relative height Eb ~ E a = Ea 
in order for the second opening to occur. The two proba- 
bilities po(t) and p a (t), obtained by numerical integration 
of eqs. Ijl Op and pip with r = 1 pN/s, are shown in fig. El 
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Fig. 2. Probabilities po and p a as functions of the time, ob- 
tained by numerical integration of eqs. 1101 and lilt with r = 1 
pN/s. 



value of r here considered (r = 10~ 2 pN/s), and at inter- 
mediate force values, q(f) lies on a line given by eq. |3"| 
with Ax — xb and AE = Eb- 

These results are in agreement with those of the above 
mentioned work of Bartolo et al. [T3|, where the authors 
found out that slope of the breaking force as a function 
of In r is equal to the position of the outermost barrier, in 
the small r regime. The same quantity has been found to 
be equal to the position of the innermost barrier at large 
r, if the relative height of the two barriers is similar, as 
in our case. It is worth to remark that in this cited work 
the single escape rate approximation has been used, i.e., 
it has been assumed that the mean escape time from each 
of the landscape barrier is constant. In the present work 
we do not use such approximation, and thus our results 
can be considered more general of those contained in • 
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Fig. 3. Plot of q as a function of /, as defined by eq. for 
different values of r. At large force, for each value of the loading 
rate r, the curves all collapse on a single scaling function which 
corresponds to the escape from the innermost barrier (xa,Ea) 
(full line). For small values of r, and in the small force regime, 
the curves also converge to a function which corresponds to the 
escape from the outermost barrier (xb,Eb) (dashed line): in 
this limit the structure of the inner part of the energy landscape 
does not affect the behavior of the function q(f). 



3 Results 

In figure the function q, as defined by eq. J2J> is plotted 
as a function of / for the energy landscape represented 
in figure ^ with the set of parameter values S given by 
eq. (|13p . and for different values of the loading rate r: the 
behavior of such function turns out to be dependent on 
both the force range and the loading rate. For all the value 
of r here considered, in the large force regime / > /„, all 
the curves collapse on a single scaling curve, which cor- 
responds to the escape from the innermost barrier Ea at 
xa, i.e., after the outermost barrier has disappeared, the 
innermost barrier becomes the only obstacle for the un- 
folding of the molecule. On the other hand, for the smallest 



The plot in fig. El indicates that the behavior of q(f), 
in the intermediate force range, strongly depends on the 
value of the pulling rate r. By using linear fits, we find 
that for all the curves with r > 1 pN/s, in the force range 
10 pN < / < f a ~ 40 pN, the slope of q(f) is equal to 
x m — xb — x a + xa = 1.6 nm. This indicates that the 
two bonds here considered, behave as a single macro-bond 
whose typical length is x m . In analogy with the single 
bond case, in the force range 10 pN < / < f a , the quantity 
q can be written as 



q(f, r) = In 



k(r) 



PfXr, 



(14) 



where we have explicitly taken into account the depen- 
dence on r, and where k(r) is the zero-force transition rate 
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of the macro-bond, which depends on the details of the en- 1 

ergy landscape. Such a quantity decreases as the loading 

rate r increases, and it turns out to scale as k(r) oc r" 1 , as 

can be seen in figure^] Using qualitative arguments, Evans 

[TU] has proposed that in a set of N identical molecular 

bonds in series, the zero force transition rate decreases as 06 

ss 

the inverse of N. In our case the inverse proportionality is 
of the zero-force transition rate on r, appears to depend § 04 
strongly on the model used here, and in particular on the 
form of eqs. (110ft and (|1 l|l used to obtain the quantity 4>(t). 
The dependence of q on the loading rate r, indicates that 



+ 





" nnn .nn» nn " 
" □□□ nDD 








t = 1 pN/s □ 




r = 10 pN/s 


¥ 


t = 10 2 pN/s - 




r = 10 s pN/s v 
ln(fc(r = l)//3xm)+0fx m - - 



I (PN) 

Fig. 4. Plot of q(f, r) + ln(r) as a function of /. The data all 
collapse on the same curve r = 1, in the force range / < f a . 



at a given force f < f a , the fraction of molecules which 
remain in the state x — before unfolding completely, in- 
creases with the loading rate r. This is also confirmed by 
a direct plot of the fraction of molecules that are in the 
ground state at the cross-over force f a , as a function of 
r, see fig. [S] In other words as r increases, the molecule 
is more and more "frozen" in its ground folded state, and 
only after the outmost barrier disappears (/ = f a ), the 
system unfolds in a way which is determined only by the 
inner barrier features. 

The characteristics of the quantity q(f, r) in the mod- 
erate force regime, i.e., the value of the slope x m — xb — 
Xa+XA and the scaling law of the zero-force transition rate 
fc(r), turn out to be fairly universal: we found the same 
results with different choices of parameter sets S (data not 
shown). This is at variance with the results discussed in 
[111113) where the slope of the rupture force as a function 
of lnr, in a given range of r, has been found to be charac- 
terized by a single length which corresponds to the relative 
position of the dominant barrier with respect to the adja- 
cent energy minimum position. In other words, our results 
suggest that the unfolding process is not dominated by a 
single escape route over a well determinate barrier, but 
is rather controlled by the interactions between the two 
barriers, which determine a different escape route. 

The analysis of the fraction of folded molecules (f>(t) 
and its related function (?(/), that we have proposed so 
far, does not provide any estimate of the intermediate en- 
ergy well E a . In the following we propose a method to ob- 



r (pN/s) 

Fig. 5. Fraction of molecules Po(fa), which "survives" in the 
state x = at the crossover force /„, as a function of r. The 
dotted line is a guide to the eye. 



tain the value of this quantity, once the fraction of (com- 
pletely or partially) folded molecules 4> and the fraction 
of molecules which are in the intermediate state p a have 
been experimentally determined, as functions of the time 
(force). Let us consider the smallest value of the pulling 
rate we have used in this paper, r = 10~ 2 pN/s. For this 
value, we have shown that 



q(f,r 



= io-) =hL («L\ 



\(3x B J 



f3E B + PfxB 



(15) 



which holds in low-to-moderate force range, see fig. 01 Us- 
ing the definition of q(f) as given by eq. J2J, we then 
obtain 



<f>(t,r = 1CT 2 ) = cxp 



_V0_ e -P(E E 



(3rx B 



-rtxs) 



(16) 



and 



fa r = Ur 2 ) = -LO e-^ EB - rtXB U(t, r = 1(T 2 ) . (17) 

The last expression is expected to hold for sufficiently 
small values of the force. On the other hand, summing 
eq. (jlUII and eq. Ijlll) . we obtain, for t < t a , 



(18) 



where fc a -*oo is given by eq. 10 • Putting together eq. (|17fl 



and eq. (|TH|) . yields 

(f>(t,r = 1(T 2 ) 



P«(t) 



^f3{E a -fx a ) 



(19) 



t=f/r 



Thus we expect that, at low-to-moderate forces where the 
equality 1)16(1 holds, the ratio of 4> to p a , for small values 
of the pulling rate, is a linear function of the force, in 
a linear-log plot of the data. This is confirmed by fig. |BJ 
where the ratio of 4>/p a , as obtained by numerical solution 
of eqs. (|TU)l and |[nj with r = 1CT 2 pN/s, is plotted as a 
function of /. From a linear fit of the relative experimental 
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Fig. 6. Ratio of folded molecule fraction (j> to the fraction p a of 
molecules in the state a function of the force /, with 

r = 10 -2 pN/s. The dashed line corresponds to the expected 
function in the small force regime, as given by eq. 1191 , 



data, it might then be possible to obtain estimates for x a 
and E a . 

The behaviour of the model here proposed appears to 
be rather robust with respect to changes in the values of 
the parameters. 

It is clear that, e.g., an increase in luq as well as a 
(slight) decrease in the barrier heights would yield a cor- 
responding decrease in the time scales. In the same way, 
a faster system time evolution could be obtained by, e.g., 
an increase of either x a or xb or a decrease in x a (see eqs. 
J7J), (JSJ and 10 which define the transition rates). 

We have shown the numerical results for a energy land- 
scape with Ea = (Eb — Ea). No qualitative change in the 
function q(f, r), with respect to that shown in figure^ has 
to be expected by changing slightly the value of one of the 
two barriers, while keeping the other one unchanged. Ob- 
viously, if one barrier becomes much larger than the other, 
the behaviour of the single barrier case is recovered. 

With our choice for the model parameters, we have 
implicitly assumed that the inner barrier is the one which 
survives longer. Let to be the time at which the inner 
barrier disappears, defined as to = Ea/ (txa), and let t max 
be the maximum time for the integration of the evolution 
equations (|10|l and Thus, with our parameter choice, 
we took tmax < to- This corresponds to stop the kinetic 
process before the last barrier disappears: in the absence of 
barriers the Kramer formalism makes no longer sense. We 
now want to discuss shortly the case where it is the outer 
barrier to survive longer, i.e. t a > to (and t a > t max > to). 
This can be obtained, e.g., by taking (EB—E a ) moderately 
greater than Ea- In this case, a slightly modified version 
of eqs. ifTUI) and ltTT|) has to be considered to take into 
account the disappearing of the barrier (xa, Ea) at t = to. 
The outcome for q(f, r) is similar to that shown in fig. 
but the slope at large force is xb — x a rather than xa- 
Still we believe that the choice Ea — (EB—E a ), and xa — 
{xb—Xo) is the most reasonable at this level of knowledge 
of the RNA tertiary structures, since each tertiary contact 
originates from the same physical mechanisms. 



4 Discussion 



The results shown in figures and 01 indicate that mea- 
suring q(f) at different values of r sheds light on differ- 
ent parts of the energy landscape. Direct information on 
the outermost barrier can only be obtained by measuring 
q(r, /) for very small values of r, while for high values of 
/, one obtains information on the innermost barrier. The 
model here presented also provides a method to obtain 
an experimental estimate of the intermediate energy min- 
imum. The measurement of such quantity is highly desir- 
able, since it determines the stability of the intermediate 
state with respect to molecule pulling. 

In the moderate force regime, and for r > 1 pN/s, 
the quantity q(f, r) does not give direct information on 
a single energy barrier, but rather indicates that in this 
force range there is a strong cooperativity in the unbind- 
ing process between the two kinetic barriers, that can be 
regarded as a single one. 

In conclusion, we have proposed a simple model for 
the unfolding of RNA molecules with tertiary structure 
under dynamic loading, where two distinct kinetic barri- 
ers with similar height, hinder the system from opening. 
This choice leads the fraction of folded molecule to depend 
on the pulling rate, and therefore accounts for experimen- 
tal results where apparently such a dependence has been 
observed [HJ ■ In order for our model to be checked, the me- 
chanical openings of a RNA molecule with tertiary con- 
tacts has to be performed with a wide range of pulling rate. 
The existence of an unique slope of q(f) vs. / in the inter- 
mediate force range, and a scaling law for the zero-force 
transition rate, as the one we find out here, would unam- 
biguously indicate that the single barrier picture is inade- 
quate to describe the kinetic process corresponding to the 
mechanical breaking of an RNA tertiary structure. Our 
results suggest that the predominance of one of the two 
single kinetic barriers on the unfolding process, cannot be 
inferred by measuring the fraction of unfolded molecules 
in a relatively large range of force. On the contrary, the ob- 
servation of this quantity over such an intermediate force 
range, suggests the existence of a complex kinetic bar- 
rier, which has a non trivial and unexpected connection 
with the single barriers. According to our model, in an 
unfolding experiment, a direct insight into the details of 
the single barriers can only be obtained either using a rel- 
atively small value of the pulling rate, or analyzing the 
unfolding kinetics in the large force regime. However the 
first possibility is limited by the need to keep the effect of 
apparatus drift under control. 
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